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Abstract 



The result of a search for heavy long-lived charged particles produced in pp colli- 
sions at i/s = 7 TeV at the LHC is described. The data sample has been collected 
using the CMS detector and corresponds to an integrated luminosity of 5.0 fb _1 . The 
inner tracking detectors are used to define a sample of events containing tracks with 
high momentum and high ionization energy loss. A second sample of events, which 
have high-momentum tracks satisfying muon identification requirements in addition 
to meeting high-ionization and long time-of-flight requirements, is analyzed inde- 
pendently. In both cases, the results are consistent with the expected background 
estimated from data. The results are used to establish cross section limits as a func- 
tion of mass within the context of models with long-lived gluinos, scalar top quarks 
and scalar taus. Cross section limits on hyper-meson particles, containing new ele- 
mentary long-lived hyper-quarks predicted by a vector-like confinement model, are 
also presented. Lower limits at 95% confidence level on the mass of gluinos (scalar 
top quarks) are found to be 1098 (737) GeV/c 2 . A limit of 928 (626) GeV/c 2 is set for a 
gluino (scalar top quark) that hadronizes into a neutral bound state before reaching 
the muon detectors. The lower mass limit for a pair produced scalar tau is found to 
be 223 GeV/c 2 . Mass limits for a hyper-kaon are placed at 484, 602, and 747 GeV/c 2 for 
hyper-p masses of 800, 1200, and 1600 GeV/c 2 , respectively. 
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1 Introduction 

Various extensions to the standard model (SM) of particle physics allow for the possibility that 
as-yet-undiscovered massive (>100 GeV/c 2 ) elementary particles could be long-lived with life- 
time greater than ~1 ns as a result of a new conserved quantum number, a kinematic constraint 
or a weak coupling (TJ^J. Such particles, where they are electrically charged, are referred to as 
Heavy Stable Charged Particles (HSCP) in this article. Because of their high mass, a significant 
fraction of the HSCPs that could be produced at the Large Hadron Collider (LHC) are expected 
to be detectable as high momentum (p) tracks with an anomalously large rate of energy loss 
through ionization (dE/ dx) and an anomalously long time-of-flight (TOF). 

Previous collider searches for HSCPs have been performed at LEP BHZI, HERA [8j, the Teva- 
tron M9til5j, and the LHC IJT6H2T) . HSCPs are expected to reach the outer muon systems of the 
collider detectors even if they are strongly interacting. In that case it is expected that a bound 
state (R-hadron) is formed in the process of hadronization ||22l - |24 ] and that the energy loss oc- 
curs primarily through low momentum transfer interactions ITT1 I25] - |27| , allowing the R-hadron 
to traverse an amount of material typical of the calorimeter of a collider experiment. However, 
the nuclear interactions experienced in matter by an R-hadron may lead to charge exchange. A 
recent study [28] of the modeling of nuclear interactions of HSCPs favours a scenario in which 
the majority of the R-hadrons containing a gluino, g (the supersymmetric partner of the gluon), 
or a bottom squark would emerge neutral in the muon detectors. Given the large uncertain- 
ties in the nuclear interactions experienced by R-hadrons, experimental strategies that do not 
rely on a muon-like behavior of HSCPs are important. Two strategies already employed are 
to search for very slow (/3 < 0.4) R-hadrons brought to rest in the detector |Tl2l[T6ll2T| , and to 
search using only inner tracker and/ or calorimeter information ||T7l - [l"9l . 

In this article we present a search for HSCPs produced in pp collisions at y/s = 7 TeV at the 
LHC, recorded using the Compact Muon Solenoid (CMS) detector II29II . The search is based on 
a data sample collected in 2011, corresponding to an integrated luminosity of 5.0 fb _1 . Events 
were collected using either of two triggers, one based on muon transverse momentum (px) 
and the other based on the missing transverse energy (E™ ss ) in the event. This event sample 
was then used for two separate selections. In the first, the HSCP candidates were defined as 
tracks reconstructed in the inner tracker detector with large dE/ dx and high pj. In the second, 
the tracks were also requested to be associated with identified muons that had a long time- 
of-flight as measured by the muon detectors. The first selection is largely insensitive to the 
uncertainties in modeling the R-hadron nuclear interactions. For both selections, the mass of 
the candidate was calculated from the measured values of p and dE/ dx. This analysis extends 
our previously published result fF7\ through the use of a larger dataset, muon TOF information, 
and track isolation requirements. This new analysis also probes additional signal models. 

2 Signal benchmarks 

The results of this search have been interpreted within the context of several theoretical models. 
Supersymmetric models [30. [33 can in some cases allow for HSCP candidates in the form 
of gluinos, scalar top quarks (stops, t\, the supersymmetric partner of the top quark), and 
scalar taus (stau, %, the supersymmetric partner of the t). We also consider a new model 
that postulates a QCD-like confinement force between new elementary particles (hyper-quarks, 
q) 11321 and allows for long-lived hyper-mesons. 

In order to study the uncertainties related to the underlying production processes, samples 
were produced with three different multiparton interaction (MPI) models [33 1: DW with CTEQ5L 
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2 Signal benchmarks 



parton distribution functions (PDF) [34] (used in IH9T ). D6T with CTEQ6L1 PDF E3 (used 
in IlZl), and Z2 with CTEQ6L1 PDF. The latter model features a harder initial-state radiation 
spectrum. The final results of this analysis are obtained with samples using the D6T MPI, which 
yield the most conservative signal selection efficiency of the three choices. 

Gluino and stop production were modelled as pair production over the particle mass range 
130-1200 GeV/c 2 using pythia v6.422 [36|. For g production, we set the squark masses to very 
high values (>7 TeV) to reflect the scenario of split supersymmetry [371 [38j. The fraction / 
of produced g hadronizing into a g-gluon state (R-gluonball) is an unknown parameter of the 
hadronization model and affects the fraction of R-hadrons that are produced as neutral parti- 
cles. As in [I7|, results were obtained for two different values of /, 0.1 and 0.5. Unless specified 
otherwise, the value / = 0.1 is assumed. The interactions of the HSCPs using the CMS appa- 
ratus and the detector response were simulated in detail with the Geant4 v9.2 f39ll40"| toolkit. 
Two scenarios for R-hadron strong interactions with matter were considered: the first follows 
the model defined in |27l |4T| . while the second is one of complete charge suppression, where 
any nuclear interaction of the R-hadron causes it to become neutral. In the second scenario, 
effectively all R-hadrons are neutral by the time they enter the muon system. 

The minimal gauge-mediated supersymmetry breaking (GMSB) model [42J was used to de- 
scribe ?i production, which can proceed either via direct pair production or via production 
of heavier supersymmetric particles that decay to one or more X\. The latter process has a 
larger cross section than direct production. Two benchmark points on the Snowmass Points 
and Slopes line 7 |43| have been considered. They correspond to N = 3 chiral SU(5) multi- 
plets added to the theory at the scales F = 100 and 160 TeV [42j respectively, and an effective 
supersymmetry-breaking scale A = 50 and 80 TeV respectively. Both points have the ratio 
of neutral Higgs field vacuum expectation values tan f> = 10, a positive sign for the Higgs- 
Higgsino mass parameter (sgn(^) = 1), and the ratio of the gravitino mass to the value it would 
have if the only supersymmetry breaking scale were that in the messenger sector, c grav = 10 4 . 
The particle mass spectrum and the decay table were produced with the program ISASUGRA 
version 7.69 [44J. The resulting X\ masses are 156 and 247GeV/c 2 , and the squark and gluino 
masses are about 1.1 and 1.7TeV/c 2 , respectively. Additional mass points in the range 100 to 
500 GeV/c 2 were obtained by varying the A parameter in the range 31 to 100 TeV and keeping 
the ratio of F to A equal to 2. In addition, direct % pair production samples were generated 
separately. 

As mentioned above, another HSCP benchmark considered in this paper is a vector-like con- 
finement model that postulates a QCD-like confinement force between new elementary par- 
ticles (hyper-quarks, q) [32l| . The hyper-quarks can be confined into SM hadron-like hyper- 
mesons such as hyper- 7T (ft), hyper-K (K), or hyper-p (p). We assume a simplified model (sim- 
ilar to that in Section 4.2 of Ref. |32|) with ft or K pair production via either the Drell-Yan 

process (K K) or via production of a resonant p (p — > K K) analogous to QCD p meson produc- 
tion [45 1 . The p mixes only with the electro weak gauge bosons, and therefore is not produced 
strongly. In this model, the ft is short-lived and decays to SM gauge bosons (e.g. ft — > W ± 7, 
ft — > W ± Z). Its production processes are not included in the simulation. However, the K is 
long-lived compared to the detector size and constitutes an HSCP candidate. In the considered 
model, the K, like the T\, would only interact via the electroweak force. The mix of resonant 
and Drell-Yan production results in different kinematics, compared to the GMSB f\ model. For 
K mass values much less than half the p mass, the HSCP receives a significant boost from the 
resonance decay [46], while near threshold the K is slow enough that the acceptance drops dra- 
matically. Parton-level events were generated with CALCHEP v2.5.4 [47J and passed to PYTHIA 
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for hadronization and simulation of the underlying event. The masses of the K, ft, and p are 
treated as free parameters in the model, affecting in particular the p width. We used a fixed ft 
mass of 600 GeV/c 2 and K masses in the range 100 to 900 GeV/c 2 , for p masses of 800, 1200 and 
1600GeV/c 2 . 

In all simulated samples, the primary collision event was overlaid with simulated minimum- 
bias events to reproduce the distribution of the number of inelastic collisions per bunch cross- 
ing (pile-up) observed in data. 

3 The CMS detector 

A detailed description of the CMS detector can be found elsewhere [29 j. The central feature 
of the CMS apparatus is a superconducting solenoid of 6 m internal diameter. Within the field 
volume are the silicon pixel and strip inner tracking detectors, the crystal electromagnetic cal- 
orimeter (ECAL), and the brass/ scintillator hadron calorimeter (HCAL). Muons are measured 
in gas-ionization detectors embedded in the steel return yoke. In addition to the barrel and 
endcap detectors, CMS has extensive forward calorimetry. CMS uses a right-handed coordi- 
nate system, with the origin at the nominal interaction point, the X axis pointing to the cen- 
ter of the LHC ring, the y axis pointing up (perpendicular to the LHC plane), and the z axis 
along the counterclockwise-beam direction. The polar angle, 9, is measured from the positive 
z axis and the azimuthal angle, (p, is measured in the x-y plane. The muons are measured in 
the pseudorapidity (tj = — lntan(0/2)) range \rj\ < 2.4, with detection planes made using 
three technologies: drift tubes (DT), cathode strip chambers (CSC), and resistive plate cham- 
bers (RPC). The DT and CSC detectors are installed in the barrel at \t] | < 1.2 and in the endcaps 
at 0.9 < \rj | < 2.4, respectively, whereas RPCs cover the range \ t] | < 1.6. The inner tracker mea- 
sures charged particle trajectories within the pseudorapidity range | rj \ < 2.5. It consists of 1440 
silicon pixel modules and 15 148 silicon strip modules. The pj resolution for tracks measured 
in the central (forward) region of the silicon tracker is 1% (2%) for pj values up to 50 GeV/c and 
degrades to 10% (20%) at pj values of 1 TeV/ c. The CMS trigger consists of a two-stage system. 
The first level (LI) of the CMS trigger system, composed of custom hardware processors, uses 
information from the calorimeters and muon detectors to select a subset of the events. The 
High Level Trigger (HLT) processor farm further decreases the event rate from around 100 kHz 
to around 300 Hz, before data storage. 

3.1 dE/dx Measurement 

The dE/dx measurement for a candidate track was performed using the charge information 
contained in the track measurements provided by the silicon strip and pixel detectors. A silicon 
strip or pixel measurement consists of a cluster of adjacent strips or pixels with a charge above 
threshold. These clusters form the basis for dE/ dx measurements in this analysis. For dE/ dx 
measurement purposes, a 'cleaning procedure' was applied to the clusters found in the silicon 
strip detectors. This selection is intended to reduce anomalous ionization contributions due to 
overlapping tracks, nuclear interactions, and hard <5-rays in the silicon strip detectors. Genuine 
single particles release charge primarily within one or two neighbouring strips. Other strips 
generally carry only a fraction (to a first approximation equal to 10~", where n is the distance 
in units of strips) of the total cluster charge from capacitive coupling and cross- talk effects [48 J. 
Measurements displaying multiple charge maxima or more than two adjacent strips containing 
comparable charge were therefore not used in the dE/ dx calculations. This cleaning procedure, 
which discards on average about 20% of the track measurements, rejects background at high 
dE/ dx without a significant impact on the signal acceptance. 
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3 The CMS detector 



As in [[17), a modified version of the Smirnov-Cramer-von Mises l49l |50| discriminant was 
used for estimating the degree of compatibility of the observed charge measurements with 
those expected for particles close to the minimum of ionization: 




Pi x Pi 



2i 
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(1) 



where / is the number of track measurements in the silicon-strip detectors, P, is the probability 
for a particle close to the minimum of ionization to produce a charge smaller or equal to that of 
the z'-th measurement for the observed path length in the detector, and the sum is over the track 
measurements ordered in terms of increasing P,-. The charge probability density functions used 
to calculate P, were obtained using reconstructed tracks with pj > 5GeV/c in events collected 
with a minimum bias trigger. As in [17], the most probable value of the particle dE/ dx was 
determined using a harmonic estimator 1^: 

-I N \ -1/2 

where N is the total number of track measurements in the pixel and silicon-strip detectors and 
c, is the charge per unit path length of the z'-th measurement. As implied above, the 1^ estimator 
was computed using both silicon strip and pixel measurements, whereas the J as estimator was 
based only on silicon strip measurements. As in [17], the mass measurement was based on the 
formula: 

2 

m 

I h = K^ + C, (3) 

where the empirical parameters K = 2.559 ± 0.001 MeV cm" 1 c 2 and C = 2.772 ± 0.001 MeV cm" 1 
were determined from data using a sample of low-momentum protons [51j. Equation ^ re- 
produces the Bethe-Bloch formula |52l with an accuracy of better than 1% in the range 0.4 < 
j6 < 0.9, which corresponds to 1.1 < (d£/dx)/(d£/dx)Mip < 4.0, where (d£/dx)Mip is the 
ionization energy loss rate of a particle at the minimum of ionization. Equation ^ implicitly 
assumes that the HSCP candidates have unit charge. For HSCPs with masses above 100 GeV/c 2 , 
the mass resolution is expected to degrade gradually with increasing mass. This effect is due 
to the deterioration of the momentum resolution and to the limit on the maximum charge that 
can be measured by the silicon strip tracker analogue-to-digital converter modules, which also 
affects the mass scale. These effects are modeled in the simulation. For an HSCP with a mass 
of 300GeV/c 2 ,the mass resolution and the most probable reconstructed mass were found to be 
16% and 280 GeV/c 2 , respectively. 

For each reconstructed track with momentum p as measured in the inner tracker, J as , and m 
were computed using Eq. p}, (0, and These values were used in the candidate selection as 
described in Sections H]ancl|5l 



3.2 TOF Measurement 

A major addition to this analysis with respect to that in Ref. Itl7l is the use of TOF informa- 
tion from the muon system. The j6 measurement for a candidate track was performed using 
time information provided by the individual DT and CSC track measurements. The DT sys- 
tem consists of four layers of muon chambers interleaved with the return yoke of the solenoid. 
The chambers in the three innermost layers contain three super-layers (SL) each, two of them 
measuring the track (p projection and the third measuring the 9 projection. The chamber in the 
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outermost layer is equipped with just two SLs that measure the track <p projection. Each SL is 
composed of four DT layers. For TOF measurement purposes, only SLs providing measure- 
ments in the (p projection were used because their time resolution is a factor of two better than 
that of the 0-projection SLs. The CSC system comprises four layers of chambers at increasing 
|z| positions. Each chamber contains six detection layers. All detection layers were used for 
TOF measurement purposes. 

Both the DTs and CSCs measure the difference (S t ) between the particle crossing time and the 
average time at which a high-momentum muon, produced at the nominal collision point in the 
triggered bunch crossing, would pass through the same portion of the detector. Measurements 
from prompt HSCPs would yield a 3t greater than zero. The (^-projection DT measurements 
within a chamber were fitted with a straight line. In order to improve the accuracy of the pa- 
rameters of the straight line for late tracks, a time shift common to all measurements within 
the chamber was introduced as a third free parameter of the fit. Having four chambers with 
eight layers measuring the track (p projections, there are up to 32 independent DT St measure- 
ments along a candidate track. Each detection layer in a CSC chamber has a nearly orthogonal 
layout of anode wires and cathode strips. The arrival time of the signals from both the anode 
wires and cathode strips measures the particle St- Having four chambers, six detection lay- 
ers per chamber, and two 5 t measurements per layer, there are up to 48 independent CSC 5 t 
measurements along a candidate track. 

A single 5 t measurement can be used to determine the track p>~ 1 via the equation: 

r x = i + ^ (4) 

where L is the flight distance and c is the speed of light. The track /3 _1 value was calculated as 
the weighted average of the /3 _1 measurements associated with the track. The weight for the 
i DT measurement is given by: 

-, = ^-£ (5) 
n Cjjj 

where n is the number of (p projection measurements found in the chamber from which the 
measurement comes and c D j- is the time resolution of DT measurements, for which the mea- 
sured value of 3 ns is used. The factor (n — 2) fn arises from computing measurement residuals 
in a plane and with respect to a straight line resulting from a fit to the same measurements. The 
weight for the i th CSC measurement is given by: 

L 2 

Wi = (6) 
°i 

where c„ the measured time resolution, is 7.0 ns for cathode strip measurements and 8.6 ns for 
anode wire measurements. 

To reduce the impact of outliers, which are mostly observed in the CSC anode wire measure- 
ment distribution, the CSC /3 _1 measurement whose difference with the track averaged fr x is 
largest was discarded if the difference was greater than three times the estimated uncertainty 
(cTg-i) in the track averaged f>~ x . The track averaged and the associated uncertainty were 
recomputed without the excluded measurement and the procedure was iterated until no fur- 
ther measurements could be discarded. A Gaussian fit to the core of the distribution of the 
measurements for the candidates passing the muon-like track pre-selection defined in Section|4] 
yielded a width of approximately 0.06, independent of the candidate pseudorapidity. 
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4 Trigger and data selection 



4 Trigger and data selection 



Events were selected using a trigger requiring a muon with high transverse momentum (pj > 
40GeV/c) with \rj\ < 2.1, or a trigger requiring large missing transverse energy (E™ ss > 150 
GeV). The latter quantity was computed online using jets reconstructed with a particle-flow 
algorithm |53) . Jet clustering was performed using the anti-Zcj algorithm Ij54j with a size pa- 
rameter of 0.5. Triggering on E™ ss allows the recovery of events with HSCPs failing muon 
identification or emerging mainly as neutral particles after traversing the calorimeters. The LI 
muon trigger accepts tracks that produce signals in the RPC detectors either within the 25 ns 
time window corresponding to the collision bunch crossing, or within the following 25 ns time 
window. This operation mode is particularly suited for detecting late tracks in the muon sys- 
tem. It was designed to cater for this analysis, and is tenable as long as collisions are separated 
by 50 ns or more, which was the case for the 2011 LHC running period. The DT and CSC LI 
triggers were used only for detecting particles produced in the collision bunch crossings. Track 
reconstruction in the muon HLT assumes particles traveling at the speed of light and produced 
within the triggered bunch crossing. However, the requirements on the quality of the muon 
segments are loose enough to allow tracks from late particles to be reconstructed with reason- 
ably high efficiency. Events with pair produced T\, and with the fastest %\ having /3 as low as 
0.6, would be selected by the muon trigger with 75% efficiency. The muon trigger efficiency 
would become less than 10% for events where the fastest X\ had f> < 0.45. 

The analysis made use of two offline selections referred to as 'tracker only' and 'tracker+TOF'. 
In the tracker-only selection, HSCP candidates were defined as individual tracks reconstructed 
in the inner tracker with large dE/dx and pj. In the tracker+TOF selection, the track was 
additionally required to be associated with an identified muon with long TOE Events selected 
online with either of the muon or E?f iss triggers were used in each of these two offline selections, 
to maximize the acceptance for HSCP signals. As described in section |6j the uncertainty in the 
signal acceptance arising from the uncertainty in the trigger efficiency is also reduced for some 
of the signals, because of the overlap of the two triggers. The tracker+TOF selection is not a 
subset of the tracker-only one because looser criteria on dE/ dx and pj can be requested in the 
former. 

For both offline selections, candidates were preselected by requiring pj (as measured in the 
inner tracker) to be greater than 45GeV/c, the relative uncertainty in the pi to be smaller than 



0.25, the track fit ^ 2 /ndf < 5, \t]\ < 1.5, and the impact parameter ^jd\ + d\ y < 0.5 cm, where 



d z and d xy are the longitudinal and transverse impact parameters with respect to the recon- 
structed primary vertex that yields the smallest track d z value. The t] requirement results from 
a search optimization based on the best discovery reach, described in section [5] Candidates 
were required to have at least two measurements in the silicon pixel detectors and at least 
eleven measurements in the inner tracking detectors before the cleaning procedure. No more 
than 20% of the inner tracker layers were allowed to be missing between the first and last mea- 
surements of the track. Candidates were required to have at least six silicon strip measurements 
passing the cleaning procedure criteria and, therefore, used for the dE/ dx and mass measure- 
ments. Candidate tracks were required to have 1^ > 3MeVcm _1 for the initial selection. For 
the tracker+TOF candidates, the additional requirements of /3 _1 > 1, where /3 was computed 
from the TOF, and a^-i < 0.07 were applied. The number of independent measurements used 
for the TOF computation was required to be greater than seven. Track candidates were required 
to be loosely isolated as measured by both the inner tracker and the calorimeters. Inner tracker 
isolation was established by considering all tracks whose direction had a distance from the can- 
didate track direction, AR = y (A<p) 2 + (A?/) 2 < 0.3. The scalar sum of the pj of these tracks, 
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with the exception of the candidate track, was required to be less than 50 (100) GeV/c for the 
tracker-only (tracker+TOF) selection. Calorimeter isolation was defined as the ratio between 
the sum of the energies measured in each ECAL and HCAL tower within a distance AR < 0.3 
from the candidate direction, and the candidate momentum. This ratio was required to be less 
than 0.3 (0.6) for the tracker-only (tracker+TOF) selection. 

Good separation between HSCPs and SM particles may be achieved by selecting candidates 
with high pj, high dE/ dx, and long TOF (in the tracker+TOF selection). These quantities are 
expected to be uncorrelated for SM particles, while a slow-moving HSCP would have high 
dE/ dx and long TOF even at high pj. Figure [T] shows the strong discrimination possible be- 
tween simulated signals and background using I as , TOF, and pj. Because of the limited number 
of available simulated QCD multi-jet events with low transverse-momentum transfers, which 
contribute to the MC distributions for SM processes, these distributions display bin-to-bin vari- 
ations in the size of the statistical errors. A disagreement was found in the tails of the I as and 
j6 _1 distributions between the data and the simulation. The I as discrepancy is understood as 
due to an increase with time of the average signal charge observed in the silicon strip detectors 
during the 2011 running period. These discrepancies in the tails have no impact on the esti- 
mated background rate since the latter is determined from data, as described in the following 
section. Signal acceptance is instead estimated from MC and studies were performed to assess 
the systematic uncertainty arising from the accuracy of the simulation model of I as and TOF. 
They are detailed in section [6] 

5 Background determination and search optimization 

The search was performed as a counting experiment in a mass window that depended on the 
HSCP mass hypothesis, M, and the model of interest. For a given M, the mass window ex- 
tended from M r eco - 2cr to 2 TeV/c 2 , where M reco is the average reconstructed mass for an HSCP 
of mass M and a is the mass resolution expected at the true HSCP mass M. The values of M reC o 
and a as a function of M were obtained from simulation. 

The candidates passing the pre-selection described in Section |4] were used for both the signal 
search and background estimate. For the tracker-only selection, signal candidates were re- 
quired to have J as arid pj greater than threshold values optimized for each model and mass 
point, as described at the end of this section. A method that exploits the non-correlation be- 
tween the pj and dE/ dx measurements for SM particles was used to estimate the background. 
The number of candidates expected to pass both the final pj and J as thresholds was estimated 
as D = BC/ A, where A is the number of candidates that fail both the J as and pj selections 
and B (C) is the number of candidates that pass only the J as (px) selection. The B (C) candidates 
were then used to form a binned probability density function in 1^ (p) for the D candidates such 
that, using the mass determination (Eq. |3]l), the full mass spectrum of the background in the 
signal region D could be predicted. It was observed that the t] distribution of the candidates at 
low dE/ dx differs from the distribution of the candidates at high dE/ dx. This effect is due to 
the typical number of measurements attached to a track, which is //-dependent and is anticor- 
related with both 7 as and 1^. The t] dependence of dE/ dx can bias the shape of the predicted 
background mass spectrum in the signal region because the p distribution is also //-dependent. 
To correct for this effect, events in the C region were weighted such that their t] distribution 
matched that in the B region. 

For the tracker+TOF selection, this method was extended to include the TOF measurement, 
assuming a lack of correlation between the TOF, pi, and dE/dx measurements. With three 
independent and uncorrelated variables, the number of background candidates in the signal 
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5 Background determination and search optimization 



CMS \fs = 7TeV 5.0 fb 1 CMS \/s = 7TeV 5.0 fb 1 
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Figure 1: Normalized distributions of pj, I as , and f> in data, simulated SM processes, and 
some of the simulated signal samples. The two plots on the left are for the tracker-only selec- 
tion. The three plots on the right are for the tracker+TOF selection. Different simulated signal 
samples are used for the left and right plots. 
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region may be estimated using six independent combinations of three out of the eight exclu- 
sive samples, each characterized by candidates passing or failing the three thresholds. These 
eight samples are analogous to the A, B, C, and D samples defined above. An additional inde- 
pendent background estimation in the signal region D may be obtained with a combination of 
four out of the eight samples. The corresponding expression is D = AGF /E 2 , where E is the 
number of candidates that fail all selections, and A, G and F are the numbers of candidates that 
pass only the I as , pj, and TOF selection, respectively. This latter estimation has the smallest 
statistical uncertainty since the four samples are such that at most one of the three thresholds 
is passed in each of them, while the other estimations have at least one suppressed population 
sample because of the requirement for two thresholds to be exceeded. For this reason the back- 
ground estimation was taken from this combination. As in the tracker-only analysis, weights 
were applied to correct the rj distribution in the regions providing the dE/ dx and TOF binned 
probability density functions that were used to model the background from SM particles in the 
signal region. The dependence of the TOF measurement on tj for genuine relativistic muons 
is due to differences in the typical number of track measurements, the accuracy of the mea- 
surements, the incident angles of particles on the detectors, and the residual magnetic field 
in the muon chamber drift volumes. The systematic uncertainty in the expected background 
in the signal region is estimated to be 10%, from the differences observed between the four 
background estimates having the smallest statistical uncertainties. The same uncertainty was 
adopted for the tracker-only selection. The statistical uncertainty of the background estima- 
tion in either the whole signal region or in a given mass range was obtained by generating 
simulated pseudoexperiments drawn from the observed distributions in the control regions. 

A 'loose' selection is defined such that there are a relatively large number of background can- 
didates in the signal region. This selection allows a cross-check on the accuracy of the back- 
ground prediction to be performed. Table [l] reports the minimum values of pj, I as , and 
that candidates must have to pass this selection, as well as the absolute yields of the back- 
ground prediction and the observed data. Figure [2] shows agreement between the observed 
and predicted mass spectra obtained using the loose selection for both the tracker-only and 
tracker+TOF candidates. The background prediction obtained from simulation using the same 
method as for data is shown in the same figure. 



Table 1: Selections used to create the 'loose' samples with large number of events and the ex- 
pected (Exp.) and observed (Obs.) event yields. The selections are defined in terms of thresh- 
olds in pj, I as , and fi^ 1 (measured from TOF). 



Selection 


p min (GeV/c) 


rmm 

i as 


o— 'Imin 


Exp. 


Obs. 


Tk-Only 


50 


0.10 




103450 ± 10350 (syst) ± 210 (stat) 


94910 


Tk+TOF 


50 


0.05 


1.05 


88010 ± 8800 (syst) ± 290 (stat) 


72079 



The final selection thresholds on pj, J as , and TOF were optimized for each signal model and 
mass by minimizing the signal cross section value for which a discovery would be achieved, 
where discovery is defined as the expected mean significance of the observed excess being 
equal to five standard deviations with at least five observed candidates. It was verified that in 
all cases the optimized thresholds also guarantee that the expected 95% confidence level (CL) 
cross section upper limit on the considered model is at most 10% larger than the minimum 
attainable. The optimized selection thresholds and the resulting signal acceptance for some 
representative signal models are reported in Tables [2] and |3j 
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5 Background determination and search optimization 



CMS Vs = 7TeV 5.0 fb" 1 
Tracker - Only 
• Observed 

A Data-based SM prediction 
~| SM prediction (MC) 

MC - Gluino (M=600 GeV/c 2 ) 




1000 

Mass (GeV/f 2 ) 



CMS Vs = 7TeV 5.0 fb - 
Tracker + TOF 

• Observed 
| A | Data-based SM prediction 
SM prediction (MC) 
MC-Stau (M=156 GeV/c*) 




1000 

Mass (GeV/c 2 ) 



Figure 2: Distribution of the candidate mass for the loose selection defined in Table [l] for the 
tracker-only (left) and tracker+TOF (right) candidates. Shown are: data, background estimate 
from data with its uncertainty (yellow band), simulated signal (green shaded histogram) and 
background prediction from MC (blue band) using the same method as for data. 



Table 2: Results of the tracker-only analysis for some representative signal mass values (in 
GeV/c 2 ): final selections in terms of minimum values of px (in GeV/c), I as , and M reco 
(in GeV/c 2 ), signal acceptance ("Acc"), number of candidates expected from SM background 
("Exp."), number of observed candidates ("Obs."), predicted theoretical cross section ("Th. 
c"), expected median cross section upper limit at 95% CL for the background-only hypothe- 
sis ("Exp. a"), and observed 95% CL cross section upper limit ("Obs. c"). All cross section 
values are expressed in pb. 



Model 


Mass 


Ft 


imin 
1 as 


M mm 

v reco 


Acc. 


Exp. 


Obs. 


Th. tr 


Exp. a 


Obs. u 


g(/ = o.i) 


300 


60 


0.400 


180 


0.16 


0.328 ± 0.040 





6.6E+01 


3.8E-03 


3.7E-03 


g(f = 0.1) 


700 


50 


0.300 


410 


0.21 


0.089 ± 0.009 


1 


2.1E-01 


3.0E-03 


4.0E-03 


g(/ = o.i) 


1100 


120 


0.225 


570 


0.15 


0.094 ± 0.010 





3.9E-03 


4.0E-03 


3.9E-03 


g(/ = 0.5) 


300 


60 


0.400 


180 


0.086 


0.328 ± 0.040 





6.6E+01 


6.9E-03 


6.8E-03 


g(/ = 0.5) 


700 


50 


0.300 


410 


0.12 


0.089 ± 0.009 


1 


2.1E-01 


5.3E-03 


7.1E-03 


g(/ = 0.5) 


1100 


120 


0.225 


570 


0.085 


0.094 ± 0.010 





3.9E-03 


7.0E-03 


6.9E-03 


g (/ = 0.1, ch. suppr.) 


300 


60 


0.400 


180 


0.020 


0.328 ± 0.040 





6.6E+01 


3.0E-02 


3.0E-02 


g (/ = 0.1, ch. suppr.) 


700 


50 


0.325 


370 


0.045 


0.092 ± 0.010 


1 


2.1E-01 


1.4E-02 


1.8E-02 


g (/ = 0.1, ch. suppr.) 


1100 


50 


0.275 


460 


0.032 


0.085 ± 0.009 


1 


3.9E-03 


1.9E-02 


2.6E-02 


tl 


200 


60 


0.400 


130 


0.14 


1.250 ± 0.160 


4 


1.3E+01 


5.8E-03 


1.1E-02 


tl 


500 


50 


0.350 


310 


0.24 


0.126 ± 0.014 





4.8E-02 


2.4E-03 


2.3E-03 


ti 


800 


50 


0.275 


450 


0.29 


0.095 ± 0.010 


1 


1.1E-03 


2.1E-03 


2.8E-03 


ti (ch. suppr.) 


200 


70 


0.400 


120 


0.021 


1.520 ±0.202 


4 


1.3E+01 


4.0E-02 


7.2E-02 


ti (ch. suppr.) 


500 


50 


0.375 


280 


0.064 


0.102 ± 0.012 





4.8E-02 


9.1E-03 


9.1E-03 


ti (ch. suppr.) 


800 


50 


0.325 


370 


0.077 


0.092 ± 0.010 


1 


1.1E-03 


8.1E-03 


1.1E-02 


GMSB Ti 


100 


65 


0.400 


20 


0.12 


6.980 ± 0.908 


7 


1.3E+00 


1.2E-02 


1.3E-02 


GMSB Ti 


494 


65 


0.350 


300 


0.64 


0.126 ± 0.014 





6.2E-05 


9.3E-04 


9.3E-04 


pair prod, fj 


100 


70 


0.400 


40 


0.11 


4.840 ± 0.608 


6 


3.8E-02 


1.2E-02 


1.5E-02 


pair prod. Ti 


308 


70 


0.400 


190 


0.39 


0.237 ± 0.030 





3.5E-04 


1.5E-03 


1.5E-03 


K(jp(800)) 


100 


70 


0.400 


10 


0.065 


4.880 ± 0.613 


6 


1.4E+00 


1.9E-02 


2.3E-02 


K (p (800)) 


500 


50 


0.350 


320 


0.61 


0.107 ± 0.012 





2.8E-04 


9.6E-04 


9.6E-04 


K(p(1200)) 


600 


50 


0.325 


370 


0.22 


0.092 ± 0.010 


1 


2.6E-03 


2.8E-03 


3.8E-03 


KQp(1200)) 


700 


50 


0.275 


440 


0.65 


0.106 ± 0.011 


1 


6.1E-05 


9.6E-04 


1.3E-03 


K (jo (1600)) 


800 


140 


0.250 


480 


0.33 


0.118 ±0.012 


1 


2.6E-04 


1.9E-03 


2.5E-03 


K (p (1600)) 


900 


135 


0.225 


530 


0.62 


0.128 ± 0.014 





1.3E-05 


9.3E-04 


9.3E-04 
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Table 3: Results of the tracker+TOF analysis for some representative signal mass values (in 
GeV/c 2 ): final selections in terms of minimum values of pj (in GeV/c), I as , and M reco 
(in GeV/c 2 ), signal acceptance ("Acc"), number of candidates expected from SM background 
("Exp."), number of observed candidates ("Obs."), predicted theoretical cross section ("Th. 
cr"), expected median cross section upper limit at 95% CL for the background-only hypothe- 
sis ("Exp. cr"), and observed 95% CL cross section upper limit ("Obs. cr"). All cross section 
values are expressed in pb. 



Model 


Mass 


M min 
PT 


rmm 
1 as 


Imin 


M mm 


Acc. 


Exp. 


Obs. 


Th. <r 


Exp. cr 


Obs. cr 


g(/ = o.i) 


300 


55 


0.175 


1.175 


180 


0.17 


0.119 ± 0.012 





6.6E+01 


3.4E-03 


3.4E-03 


g(/ = o.i) 


700 


110 


0.050 


1.125 


430 


0.19 


0.113 ± 0.015 





2.1E-01 


3.0E-03 


3.0E-03 


g(/ = o.i) 


1100 


110 


0.025 


1.075 


620 


0.13 


0.111 ± 0.033 





3.9E-03 


4.6E-03 


4.6E-03 


g(/ = 0.5) 


300 


55 


0.175 


1.175 


180 


0.094 


0.119 ± 0.012 





6.6E+01 


6.3E-03 


6.2E-03 


g(/ = 0.5) 


700 


110 


0.050 


1.125 


430 


0.11 


0.113 ± 0.015 





2.1E-01 


5.4E-03 


5.3E-03 


g(/ = 0.5) 


1100 


110 


0.025 


1.075 


620 


0.072 


0.111 ± 0.033 





3.9E-03 


8.2E-03 


8.2E-03 


ti 


200 


50 


0.200 


1.200 


130 


0.15 


0.109± 0.011 





1.3E+01 


3.9E-03 


3.8E-03 


ti 


500 


60 


0.075 


1.150 


330 


0.25 


0.125 ± 0.013 





4.8E-02 


2.4E-03 


2.4E-03 


ti 


800 


105 


0.025 


1.125 


490 


0.26 


0.096 ± 0.019 





1.1E-03 


2.2E-03 


2.2E-03 


GMSB fi 


100 


50 


0.300 


1.275 


30 


0.20 


0.093 ± 0.011 





1.3E+00 


2.9E-03 


2.9E-03 


GMSB X\ 


494 


55 


0.025 


1.175 


320 


0.78 


0.113 ± 0.014 


1 


6.2E-05 


7.8E-04 


1.1E-03 


pair prod. Tj 


100 


50 


0.250 


1.275 


50 


0.19 


0.109 ± 0.012 





3.8E-02 


3.0E-03 


2.9E-03 


pair prod. X\ 


308 


65 


0.125 


1.200 


190 


0.55 


0.105 ± 0.011 





3.5E-04 


1.1E-03 


1.1E-03 


K(p(800)) 


100 


50 


0.300 


1.275 


20 


0.11 


0.095 ± 0.011 





1.4E+00 


5.3E-03 


5.2E-03 


K(p(800)) 


500 


60 


0.075 


1.150 


330 


0.68 


0.125 ± 0.013 





2.8E-04 


8.6E-04 


8.5E-04 


K(p(1200)) 


600 


70 


0.025 


1.150 


380 


0.22 


0.107 ± 0.015 





2.6E-03 


2.6E-03 


2.6E-03 


K(p(1200)) 


700 


110 


0.050 


1.125 


450 


0.66 


0.087 ± 0.013 





6.1E-05 


9.0E-04 


9.0E-04 


K (p (1600)) 


800 


50 


0.050 


1.100 


500 


0.33 


0.119 ± 0.021 





2.6E-04 


1.8E-03 


1.8E-03 


K (fi (1600)) 


900 


85 


0.075 


1.075 


550 


0.61 


0.123 ± 0.022 





1.3E-05 


9.3E-04 


9.1E-04 



6 Results 

After comparing data in the signal region with the expected background for all optimized se- 
lections, no statistically significant excess was observed. Tables [2] and [3] report results for some 
representative selections. The largest excess has a significance of 1.75 one-sided Gaussian stan- 
dard deviations and was found with the selection optimized for a ti with a mass of 200 GeV/c 2 . 
Only one of the three highest mass candidates passing the tracker-only loose selection (Fig. [2j 
passes one of the final selections. This candidate is also associated with an identified muon and 
has j6 _1 = 1.03, which is well below any threshold used in the tracker+TOF final selections. 

The observed data sample was used to calculate upper limits on the HSCP production cross 
section for each considered model and mass point. The cross section upper limits at 95% 
CL were obtained using a CL S approach [55| with a one-sided profile likelihood test statistic 
whose p-values were evaluated by generating pseudoexperiments using a frequentist prescrip- 
tion [56 1 . A log-normal probability density function [50| was used for the nuisance parameter 
measurements, which are the integrated luminosity, the signal acceptance, and the expected 
background yield in the signal region. When generating pseudoexperiments for the limit cal- 
culation, each nuisance parameter was drawn from the corresponding probability density func- 
tion with a central value equal to the best fit value to data under the signal+background hy- 
pothesis. All systematic uncertainties are summarized in Table |4] and are incorporated in the 
limits quoted below. 

Simulation was used to determine the signal acceptance. A number of studies were undertaken 
to estimate the degree to which the simulation correctly models the detector response to HSCPs 
and to assess an uncertainty in the signal acceptance. 

The uncertainty in the trigger efficiency derived from simulation was evaluated separately for 
the E™ ss and single-muon triggers. The uncertainty in the E™ 1SS trigger efficiency was domi- 
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6 Results 



nated by the uncertainty in the jet-energy scale l|57l , which was less than 3% across the energy 
range. For the charge-suppression models, where the E™ 1SS trigger is most relevant, varying the 
jet-energy scale and jet-energy resolution within their uncertainty resulted in a relative change 
of the trigger efficiency by no more than 5%. For all the other models, which made use of 
overlap with the single-muon trigger, the relative change of the overall trigger efficiency was 
found to vary by no more than 2%. For the single-muon trigger, a relative disagreement of up 
to 5% was observed between the efficiency estimated in data and MC at all energies Il58l . In 
addition, for this specific analysis, a further uncertainty arises from the imperfect simulation of 
the synchronization of the muon trigger electronics. The accuracy of the synchronization was 
evaluated from data separately for each muon subdetector. This effect was found to yield less 
than 2% relative uncertainty on the overall trigger efficiency for all considered signals. On the 
basis of these numbers, an uncertainty of 5% on the overall trigger efficiency was assumed for 
all models. 

The accuracy of the dE/ dx model used in simulation was studied using low-momentum pro- 
tons and kaons. The simulation was found to underestimate both the 1^ and I as scales by less 
than 5%. The I as resolution was in contrast found to be overestimated by a constant value of 
0.08 in the region around the thresholds adopted in the analysis. After corrections for these dis- 
crepancies were applied to simulation, only 20% of the signal models displayed an efficiency 
decrease, with the maximum relative reduction being smaller than 2%. The efficiency for all 
other models increased by up to 10% relative to the uncorrected MC result. Based on these 
results, the efficiency determined from the simulation was not corrected, but was assigned an 
associated uncertainty of 2%. 

The accuracy of the TOF model implemented in the simulation was studied using cosmic ray 
muons and muons produced directly in collisions. In the region around the thresholds 
adopted in the analysis, the simulation was found to overestimate by a constant value 
of 0.003 and 0.02 for DT and CSC, respectively. The resolution of the measurement of was 
found to be well modelled in simulation. After corrections for these discrepancies were applied, 
the signal efficiency was found to decrease by no more than 2% for all considered models and 
mass points. This maximum change of 2% was adopted as the uncertainty associated with the 
TOF measurement. 

The uncertainty on the track momentum scale was modelled by varying the track pj as a func- 
tion of the track (p and rj values such that |59|: 

1 1 

-r = — + 8K T (q,4>,y), (7) 



S KT (q,<p,T]) = a + brj 1 + qdsir\{(p-<po), (8) 

where q is the track charge sign (q = ±1) and the function Sk t controls the shift in the track 
momentum scale. This function has four free parameters, a, b, d, and (po, whose values were 
obtained [59 1 by minimizing the difference between the invariant mass distributions of Z — > 
candidates in data and simulation. The obtained values are a = 0.236 TeV _1 c, b = 
— 0.135 TeV _1 c, d = 0.282 TeV _1 c, and cpo = 1.337 rad. The phi dependence is believed to be 
due to imperfect inner tracker alignment. The expected shift in inverse pj for tracks of higher 
momenta measured in the inner tracker are found [59 1 to be compatible with those provided by 
equations [7] and [HJ The difference between the signal acceptance with the nominal and shifted 
pj was taken as the uncertainty and was found to be smaller than 4%. 
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Table 4: Sources of systematic uncertainties and corresponding relative uncertainties. 



Source of systematic uncert. 


Relative uncert. (%) 


i^Ji. f-yi. LU.1 U.V_V_V_- ys LUX LV_\_- • 




- Trigger efficiency 


5 


- Track momentum scale 


<4 


- Ionization energy loss 


2 


- Time-of-flight 


2 


- Track reconstruction eff. 


< 2 


- Muon reconstruction eff. 


< 2 


- Pile-up 


< 0.5 


Total uncert. in signal acc. 


7 


Expected background 


10 


Integrated luminosity 


2.2 



The uncertainties in the efficiencies for reconstructing muons [58], and for reconstructing tracks 
in the inner tracker [60J were also considered and established to be less than 2% in each case. 

The impact of the uncertainty in the mean rate of additional interactions in each bunch crossing 
was studied and found to be negligible compared to the statistical precision (0.5%) allowed by 
the size of the simulated signal samples. 

Two theoretical uncertainties affecting the signal acceptance were studied: the uncertainty in 
the model of hadronization and nuclear interactions, and the uncertainty due to the MPI tune. 
The hadronization and nuclear-interaction model is discussed in Section |2j Results are ob- 
tained for two very different nuclear interaction models and for two different g hadronization 
schemes. With regard to the MPI tune, tune Z2 uses a pj-ordered model, which appears to 
generate significantly more initial-state radiation than the Q 2 -ordered D6T model. For some 
models a significant increase in the trigger efficiency and in the reconstruction efficiency is 
found and the observed limits become more stringent. The most conservative set of limits, 
resulting from the Q 2 -ordered D6T model, are those reported. 

An uncertainty of 2.2% is estimated [61 j for the absolute value of the integrated luminosity. 
The uncertainty in the expected background was discussed in Section [5] and is estimated to be 
of the order of 10%. This uncertainty has very little impact on the results, because of the small 
numbers of expected events for most mass points. 

The 95% CL cross section upper limit curves obtained with both the tracker-only and the 
tracker+TOF selection are shown in Figs. [3] and |3J along with the theoretical predictions for 
the chosen signal models. The ratio of observed to expected 95% CL upper limits on the cross 
section is reported in Fig.|5]for the different combinations of models and scenarios considered. 
Numerical values for the predicted theoretical cross section, and expected and observed cross 
section upper limit at 95% CL are reported in Tables [2] and [3] for some representative signal 
models. For ti and g pair production, theoretical cross sections were computed at next-to- 
leading order (NLO) plus next-to-leading-logarithmic (NLL) accuracy [162TI661 using PROSPINO 
v2.1 [67]. The uncertainty in these theoretical cross section values varies between 10% to 25% 
and is shown in Fig.|3]as a band around the central value. The cross sections for the models with 
Ti production were calculated at NLO with PROSPINO v2.1. The uncertainty in the theoretical 
cross section was estimated to be 5% to 14% for the GMSB model and 3% to 7% for direct % pair 
production, depending on the mass. In all cases the sources of uncertainty include renormaliza- 
tion and factorization scales, oc s , and the parton distribution functions. The cross sections for K 
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production used in this paper are computed at leading order only. The theoretical uncertainty 
was not evaluated because of the lack of corresponding theoretical NLO calculations. For a 

fixed p mass, the K K cross section is a combination of a p resonance and Drell-Yan production. 
When the K mass is much smaller than half the p mass, Drell-Yan production dominates. As 
the K mass increases, resonance production becomes increasingly important, and dominates as 

the kinematic limit for K K pair production is approached. For K mass greater than half the p 
mass, resonance production turns off, resulting in a steep drop in the total cross section (shown 

by the nearly vertical line in Fig. ij. In addition, near the kinematic limit the p — > K K process 
produces very low velocity K particles. The signal acceptance therefore decreases dramatically 
until the resonance production turns off, at which point the acceptance increases again. This 
results in a spike in the cross section limit near the kinematic limit. 

From the intersection of the cross section limit curve and the central value of the theoretical 
cross section band, a 95% CL lower limit of 1098 (1046) GeV/c 2 on the mass of pair produced g, 
hadronizing into stable R-gluonballs with 0.1 (0.5) probability, is determined with the tracker- 
only selection. The tracker+TOF selection gives a lower limit of 1082 (1030) GeV/c 2 for the same 
signal model. The analogous limit on the ti mass is 714 GeV/ c 2 with the tracker-only selection 
and 737 GeV/ 'c 2 with the tracker+TOF selection. The charge suppression scenario discussed 
above yields a g mass limit of 928GeV/c 2 for / = 0.1 and 626GeV/c 2 for the tj. The limits 
on GMSB and pair produced X\ are calculated at 314 and 223GeV/c 2 , respectively, with the 
tracker+TOF selection. The mass limits on K are established at 484, 602 and 747 GeV/ 'c 2 for p 
masses of 800, 1200 and 1600 GeV/c 2 , respectively, with the tracker+TOF selection. 
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Figure 3: Predicted theoretical cross section and observed 95% CL upper limits on the cross 
section for the different signal models considered: production of t\, g, and T\; different fractions 
/ of R-gluonball states produced after hadronization; standard and charge suppression (ch. 
suppr.) scenario. Left: tracker-only selection. Right: tracker+TOF. The uncertainties in the 
theoretical cross section are shown as bands around the central values. 
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Figure 4: Predicted theoretical cross section and observed 95% CL upper limits on the cross 
section for the K models with different p mass values. Left: tracker-only selection. Right: 
tracker+TOF. 

7 Summary 

The CMS detector has been used to search for highly ionizing, high-pj and long time-of-flight 
massive particles. Two inclusive searches have been conducted: one that uses highly ionizing 
tracks reconstructed in the inner tracker, and a second requiring that these tracks also be iden- 
tified in the CMS muon system and have long time-of-flight. The former is model-independent 
in that it is insensitive to the details of R-hadron nuclear interactions. In each case, the observed 
number of candidates is consistent with the expected background. Upper limits on production 
cross section and lower limits on masses of stable, weakly- and strongly-interacting particles 
have been established for a variety of models. They range from 223 GeV/ c 2 for a pair produced 
scalar tau to 1098GeV/c 2 for a pair-produced gluino. These limits are the most restrictive to 
date. 

8 Acknowledgements 

We congratulate our colleagues in the CERN accelerator departments for the excellent perfor- 
mance of the LHC machine. We thank the technical and administrative staff at CERN and other 
CMS institutes, and acknowledge support from: FMSR (Austria); FNRS and FWO (Belgium); 
CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC 
(China); COLCIENCIAS (Colombia); MSES (Croatia); RPF (Cyprus); MoER, SF0690030s09 and 
ERDF (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); 
BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NKTH (Hungary); DAE and DST 
(India); IPM (Iran); SFI (Ireland); INFN (Italy); NRF and WCU (Korea); LAS (Lithuania); CIN- 
VESTAV, CONACYT, SEP, and UASLP-FAI (Mexico); MSI (New Zealand); PAEC (Pakistan); 
MSHE and NSC (Poland); FCT (Portugal); JINR (Armenia, Belarus, Georgia, Ukraine, Uzbek- 
istan); MON, RosAtom, RAS and RFBR (Russia); MSTD (Serbia); MICINN and CPAN (Spain); 
Swiss Funding Agencies (Switzerland); NSC (Taipei); TUBITAK and TAEK (Turkey); STFC 



16 



8 Acknowledgements 



CMS t^ = 7TeV 5.0 ft) -1 



CMS ^s = 7TeV 5.0 fb _ 



Tracker - Only 
Expected ± 1o 



gluino; 10% gg 



gluino; 50% gg 



gluino; 10% gg; ch. suppr. 



stop 



stop;ch. suppr. 



GMSB stau 



Pair Prod, stau 



Hyper-K,p = 0.8TeV 



Hyper-K,p = 1.2TeV 



Hyper-K,p = 1.6TeV 



500 



1000 

Mass (GeV/c 2 ) 



T3 

U 
<D 
Q. 
X 
UJ 

o 



tu 
DC 

in 
I 



2 

1- 



2 

1- 



2- 
1- 



Tracker + TOF 


gluino; 10% gg 




Expected ± 1o * ** 






— Expected ± 2o 


gluino; 50% gg 






_ 'Observed 







stop 



GMSB stau 



Pair Prod, stau 



Hyper-K,p = 0.8TeV _: 



Hyper-K,p = 1.2TeV - 



Hyper-K,p = 1.6TeV 



500 



1000 

Mass (GeV/c 2 ) 



Figure 5: Ratio of observed 95% CL upper limits to expected median limits for the background- 
only hypothesis. The green (dark) and yellow (light) bands indicate the ranges that are ex- 
pected to contain 68% and 95% of all observed excursions from the expected median, respec- 
tively. Ratios are presented for the different signal models considered: production of g, t\, X\, 
and K; different fractions, /, of R-gluonball states produced after hadronization; standard or 
charge suppression (ch. suppr.) scenario. Left: tracker-only selection. Right: tracker+TOF. 
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